Dynamical tunneling in mushroom billiards 
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We study the fundamental question of dynamical tunneling in generic two-dimensional Hamilto- 
nian systems by considering regular-to-chaotic tunneling rates. Experimentally, we use microwave 
spectra to investigate a mushroom billiard with adjustable foot height. Numerically, we obtain 
tunneling rates from high precision eigenvalues using the improved method of particular solutions. 
Analytically, a prediction is given by extending an approach using a fictitious integrable system to 
billiards. In contrast to previous approaches for billiards, we find agreement with experimental and 
numerical data without any free parameter. 

PACS numbers: 05.45.Mt, 03.65.Sq, 03.65.Xp 
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Typical Hamiltonian systems have a mixed phase space 
in which regular and chaotic motion coexist. While clas- 
sically these regions are separated, quantum mechani- 
cally they are coupled by tunneling. This process has 
been called "dynamical tunneling" pL,] as it occurs across 
a dynamically generated barrier in phase space. Tun- 
neling has been studied between symmetry related reg- 
ular regions (chaos- assisted tunneling) [2i,13, 0, S 0, 21 
and from a single regular region to the chaotic sea H 
fiol . fill . fl2| . In contrast to the well understood ID tun- 
neling through a barrier, the quantitative prediction of 
dynamical tunneling is a major challenge. Results have 
been found for specific systems or system classes only, 
e.g. recently for 2D quantum maps with an approach us- 
ing a fictitious integrable system [l2|. However, a precise 
knowledge of tunneling rates is of great importance. Re- 
cent examples are spectral statistics in systems with a 
mixed phase space |l3|, eigenstates affected by flooding 
of regular islands [IJ ^ and emission properties of optical 
micro-cavities [l5[. 

Billiards are an important class of Hamiltonian sys- 
tems. Classically, a point particle moves along straight 
lines inside a domain with elastic reflections at its bound- 
ary. Quantum-mechanical approaches for dynamical tun- 
neling rates have so far escaped a full quantitative pre- 
diction as they required fitting by a factor of 6 for the 
annular billiard [3] and by a factor of 100 (see below) for 
the mushroom billiard [16]. 

In this paper we present a combined experimental, the- 
oretical, and numerical investigation of dynamical tun- 
neling rates in mushroom billiards [17] , which are of great 
current interest [I^, [l^, [l^, [l^, [20] due to their sharply 
divided phase space. Experiments are performed using a 
microwave cavity. Extending the approach using a ficti- 
tious integrable system [l2| to billiards, we find quanti- 
tative agreement in the experimentally accessible regime. 



see Fig. (U without a free parameter. In addition, numer- 
ical computations verify the predictions over 18 orders 
of magnitude with errors typically smaller than a factor 
of 2, see Fig. SI The theoretical approach thus provides 
unprecedented agreement for tunneling rates in billiards. 

We consider the desymmetrized mushroom billiard, i.e. 
the 2D autonomous system i^(p,q) = p^/2M + ^(q) 
shown in Fig. [2)3, characterized by the radius of the quar- 
ter circle the foot width a and the foot height I. The 
potential is zero inside the domain and infinite outside. 
Classically one has regular and chaotic dynamics as vi- 
sualized by the phase-space portrait in Fig. [2]i. Quan- 
tum mechanically the billiard is described by the time- 
independent Schrodinger equation — A?/^(q) = Etlj{q) 
with Dirichlet boundary conditions at dft, using the nat- 
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FIG. 1: (color online) Dynamical tunneling rates from the reg- 
ular region to the chaotic sea for quantum numbers n < 5 vs 
kR for a mushroom billiard {a/R = 10/19): theoretical pre- 
dictions (connected by solid lines) from Eq. (|8]), experimental 
results (crosses), and numerical data (dots). The dashed line 
denotes the lower limit of the experimental resolution. 
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ural units 2M = h = 1. The eigenstates can be classified 
as being either mainly regular or mainly chaotic, depend- 
ing on the phase space region on which they concentrate. 

The initial tunneling decay of a purely (unperturbed) 
regular state, is, according to Fermi's golden rule, de- 
scribed by a rate 



7 = 2^(|^|2)pch, 



(1) 



where 



is the averaged squared matrix element be- 



tween this regular state and the chaotic states of sim- 
ilar energy. According to Weyl's formula the density 
of chaotic states is pch ^ ^ch/47r, where Ach = la -\- 
[R^ arcsin (a/ R)^a\/ W' — a^]/2 is the area of the billiard 
times the fraction of the chaotic phase-space volume. 

Under variation of the foot height / (and for not too 
large pch 13 ) avoided crossings of this regular state with 
different chaotic states are observable. The splittings 
= 2\v\ determine individual matrix elements for dif- 
ferent billiards with varying pch- We therefore extend the 
average in Eq. ^ also over pch relying on the assump- 
tion that the regular-to-chaotic tunneling rate 7 is a local 
property of the hat region, where the regular phase-space 
component resides, and thus does not depend on the foot 
height. This leads to 



7=(|Ai?|2Aeh/8), 



(2) 



where A^h in the experiment varies from 0.5 to 1.2 when 
increasing / from to 25.7 cm. 

Fig. [2^ shows the mushroom billiard used in the mi- 
crowave experiment. Spectra have been taken as a func- 
tion of the foot height / of the mushroom in the frequency 
regime 3.0 to 13.5 GHz, corresponding to values of kR 
between 11.9 and 53.8. Fig. [3] shows part of the obtained 
spectra in a small kR window. As the energy of the regu- 
lar states of the quarter circle do not depend on the foot 
height, they appear as straight horizontal lines, whereas 
the chaotic states are shifted to lower energies with in- 
creasing foot height reflecting the increasing density 
Pch of chaotic states. 

For each of the regular states, see Eq. with ra- 
dial quantum numbers n between 1 and 5, and azimuthal 
quantum numbers m even between 8 and 32 all accessi- 
ble splittings at avoided crossing A/c have been deter- 
mined by means of a hyperbola fit (Fig.[3|). From this we 
get the energy splittings = 2kAk of the correspond- 
ing quantum system and by averaging over all avoided 
crossings, Eq. ([2]), deduce the tunneling rates jmn from 
regular states (m, n) to the chaotic sea. Apart from the 
results for n = 1 they are in very good agreement with 
the theoretical prediction, Eq. ([8|), derived below. 

The experimental resolution of avoided crossings is 
limited by the line widths of the resonances caused by 
wall absorption and antenna coupling. In the studied 
frequency regime the line widths were about Au^j = 
0.01 GHz, corresponding to a Ak^R ~ 0.004. From the 
hyperbola fit of the avoided crossings all splittings AkR 
larger than O.lAk^R could still be resolved, correspond- 
ing to tunneling rates 7 between 0.001 and 0.024, see 



Fig. [T] (dashed line). Another complication is caused_W 
the antenna giving rise to an additional splitting |13( , 
which is proportional to the product of the involved 
wave functions |'0(qa)| at the antenna position q^. For 
the rightmost three data points for n = 1, see Fig. (U 
l'^reg(^a)l particularly large, which is probably the ex- 
planation for the deviations between experiment and the- 
ory observed in these cases. 

Now we derive a formula for tunneling rates in bil- 
liards. This faces a general problem: The matrix element 
V, appearing in Eq. ([1]), cannot be calculated from the 
nominally regular and chaotic eigenstates of i^, as they 
have small, but still too large, admixtures of the other 
type of states. Instead, we determine v by introducing 
a fictitious regular billiard system H^eg with purely reg- 
ular eigenstates, extending an approach for ID quantum 



maps 



i2i. 



i^reg has to be chosen such that its classical 



dynamics resembles the classical motion corresponding to 
H within the regular region as closely as possible. The 
eigenstates ?/^reg of i^reg are localized in the regular region 
and continue to decay into the chaotic sea. With states 
i/jch living in the chaotic region of phase space, we use in 
analogy to Ref. [12] for the coupling matrix element 



Jn 



reg)i^reg{x,y)dxdy. (3) 






FIG. 2: (a) Experimental desymmetrized mushroom billiard 
with radius = 19 cm, foot width a = 10 cm and foot height 
/ = . . . 25.7 cm. The antenna is located 4 cm below the top 
and has a distance of 0.75 cm from the vertical wall, (b) 
Schematic picture showing the coordinate systems used in 
the theoretical derivation, (c) Auxiliary billiard H^^. (d) 
Phase-space portrait at the quarter circle boundary (relative 
tangential momentum p vs arclength s) showing regular and 
chaotic regions with illustrations of trajectories. 
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FIG. 3: (color online) Part of the evaluated experimental res- 
onance spectra of the mushroom microwave billiard vs foot 
height /. The horizontal line corresponds to the eigenfre- 
quency of iprei^ ^ the crosses mark the extracted resonances. 
The inset shows a magnification of one avoided crossing, in- 
cluding a hyperbola fit (solid line), and the numerically ob- 
tained regular and chaotic states involved in the crossing. 



FIG. 4: (color online) Tunneling rates from regular states 
with quantum numbers n = 1, 2, and 3 vs kR for a/R — 0.5 
comparing prediction Eq. JS)) (connected by solid lines) and 
numerical data (dots). The insets show the regular eigenfunc- 
tions ^rei^{x,y) and ipref{x,y). 



Eqs. ([T]) and (|3]) define our approach for determining dy- 
namical tunneling rates in billiards. Note, that it re- 
quires the determination of reasonably good i^reg and 
tpreg^ which for a general billiard is a difficult task. 

We will now apply this approach to the desymmetrized 
mushroom billiard, see Fig. [2)3. We set R = 1 in the fol- 
lowing analysis. A natural choice for the regular system 
Hreg is the quarter-circle billiard with its eigenstates 



ij^^{r, (p) NmnJm {jmnr) sm{m(p) , 



(4) 



in polar coordinates (r, (p) . They are characterized by the 
radial (n = 1,2,...) and the azimuthal (m = 2, 4, ... ) 
quantum numbers. Here Jm denotes the mth Bessel func- 

/S/tt/ Jm-l{jmn) 



tion, jmn the nth root of J^, Nmn 
the normalization, and E^.n = Jmn eigenenergy. 



Evaluating Eq. ([3|) leads for y < to the undefined 
product of H — H^eg = — oo and i/j^^ = 0. We therefore 
introduce the auxiliary billiard i^J^, see Fig. [2t, with a 

-- W :^ E. We 
where ap- 



large but finite potential V{x^y < 0) = 
evaluate Eq. (|3]) in the limit W ^ oo, 
proaches i^reg, leading to 



V = lim 

w 



im / dx dyip^i,{x,y){-W)'tp^lw{x,y) 
dx V'ch(x, y = 0) dyi;Z^{x, y = 0), (5) 



/ 

^0 



where we performed the ^/-integration on ilj'^^y^{x^y) = 

i^7^g,wi^^y = 0)exp(V^ - EmnV), following from the 
Schrodinger equation for y < and the continuity at 
y = 0. Furthermore we used dyilj^^y^{x,y = 0) = 

— Emni^^g^wi^^ y — 0)' following from the con- 
tinuity of the derivative at = 0, and we replaced 



lim^^^oo vWex.p{\/W — EmnV) by a Dirac delta func- 
tion. Below we use dyilj^^{x,0) = mNmnJm{jmnx)/x. 



For the chaotic states tpchix^y) we employ a random 
wave description which has recently been extended 
to systems with a mixed phase space [22|. While this 
describes the behavior inside the billiard accurately, it 
would not include the effect of the boundary, e.g. near 
the corner, where the main contribution of the integral 
in Eq. ([5]) arises. W e extend a boundary-adapted ran- 
dom wave model [23] to the case of a corner with angle 
37r/2 using basis states with Dirichlet boundary condi- 
tions [2I, 



3A, 



ch 



^c^J^ (V^p) sin (^i 

s=l ^ 



(6) 



where the polar coordinates (p, 1^) at the corner are re- 
lated to (x^y) hy X = a -\- pcos{'3) and y = psm{'d) (see 
Fig. [2)3). The coefficients Cg of this ensemble are inde- 
pendent Gaussian random variables with (cg) = and 
{csCt) = Ss^t' The normalization is chosen such that 
(iV^chlp, = iMch holds far away from the corner. 
Note, that (i) we do not require these chaotic states 
to decay into the regular island, as Eq. ([5|) is an inte- 
gral along a line of the billiard where the phase space is 
fully chaotic, and that (ii) near the boundary, but away 
from the corner, one recovers the behavior 1 — Jo{2k\x\) 
[23I , [25I . Inserting Eq. (|6|) for 1? = tt into Eq. ([5]) one can 
determine the averaged squared matrix element, 
and with Eq. ([1]) one gets 



mn / _j 



a 

/(\x 
— Jm (jmnX) J ^ {jmn [a-x] ) 



(7) 

where the sum over s excludes all multiples of 3, which 
is indicated by the prime. The remaining integral can be 
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solved analytically, leading to 

8 Jm+^UrnnCl)'^ 

for the tunneling rates from any regular state V^J^g in the 
mushroom billiard. The sum has its dominant contribu- 
tion for s = 1 and using 5 < 2 is sufficiently accurate. 

It is worth to remark that a very plausible estimate 
of the tunneling rate is given by the averaged square of 
the regular wave function on a circle with radius a, i.e. 
the boundary to the fully chaotic phase space, yielding 
7mn = ^mnJmUmna)'^/^' Surprisingly, it is just about 
a factor of 2 larger for the parameters we studied. In 
Ref. [l6| a related quantity is proposed, given by the 
integral of the squared regular wave function over the 
quarter circle with radius a. This quantity, however, is 
too small by a factor of order 100 for the parameters 
under consideration. 

The eigenvalues and eigenfunctions of the mush- 
room billiard are determined by numerically solving the 
Schrodinger equation. Because of its superior compu- 
tational efficiency we have chosen to use the improved 
method of particular solutions [l6", ^26] allowing a deter- 
mination of the energies E with a relative error ^ 10~^^. 
Analyzing avoided crossings of a given regular state with 
typically 30 chaotic states we deduce from Eq. ([2]) the 
tunneling rate. Note, that some pairs of regular states 
are very close in energy, e.g. £^20,1 — ^16, 2 ^ 10~^, such 
that their avoided crossings with a chaotic state overlap. 



making a numerical determination of the smaller tunnel- 
ing rate unfeasible within the presented approach. Fig. [H 
shows the tunneling rates jrnn for fixed radial quantum 
number n = 1,2,3 and increasing azimuthal quantum 
number m, comparing the theoretical prediction, Eq. ([8]), 
with numerical results. We find excellent agreement for 
tunneling rates jrnn over 18 orders of magnitude. 

In the experiment the unavoidable coupling to the en- 
vironment dominates for small tunneling rates, but Fig.[T] 
shows, that in the microwave experiment the coupling by 
the antenna is negligible over three orders of magnitude. 
This is a promising aspect for future experimental studies 
of more complex systems, in particular when numerical 
and theoretical results are not available. 

In summary, we have presented an experimental, nu- 
merical, and theoretical investigation of tunneling rates 
in the mushroom billiard. We find agreement without 
any free parameter, which is unprecedented for billiards. 
This success of the approach using a fictitious integrable 
system gives confidence that in the future it can be ap- 
plied to generic billiards, where the determination of a 
suitable H^eg is more challenging. 
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